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ABSTRACT 


General  harmonic  solutions  have  been  found  for  the 
camber  problem  of  a  ducted  propeller  with  finite  blade 
number  in  uniform  motion  at  zero  incidence.  The  fluid  is 
assumed  inviscid  and  incompressible.  Principal  attention 
is  focused  on  the  zeroth  harmonic  based  on  a  successive 
iteration  of  the  two-dimensional  or  sectional  Glauert 
coefficients.  In  matrix  form,  the  speed  of  convergence  and 
identification  of  the  limit  are  studied.  Appropriate 
tables  for  typical  shrouds  are  provided  together  with  an 
illustrative  examplu.  Both  the  iteration  procedure  eind 
direct  simultaneous  inversion  yield  solutions  for  the  higher 
harmonics . 


iii 


TABLE  OF  CONTENTS 


INTRODUCTION 

CHAPTER  ONE  -  BASIC  FORMULATION 

1.1  Mathematical  Modex  and  Background 

1.2  Reduction  of  the  Governing  Equation 

1.3  The  Zeroth  Harmonic  Equation 

1.4  The  Higher  Harmonic  Equations 

CHAPTER  TWO  -  DETERMINATION  OF  THE  NET  AVERAGE 
DUCT  LOAD 

2.1  General  Considerations 

2.2  Expansion  of  the  Kernel  Function 

2.3  Iterative  Solution 

2.4  Convergence  of  the  Iteration 

2.5  Example 

2.6  Direct.  Inversion 

CHAPTER  THREE  -  HIGHER  HARMONIC  SOLUTIONS 

3.1  Methods  of  Approach 

3.2  Solution  by  Reduction  to  Zeroth 
Harmonic  Form 

3.3  Direct,  Simultaneous  Inversion 

3.4  comparison  of  Methods 

CONCLUSIONS 


REFERENCES 


NOMENCLATURE 


*m0  *  *^mv 


Olauert  coefficients  for  A  «  m  >  0  ; 

see  Eqs.  (3*3)  ™ 


coliimn  matrix  of  Glauert  coefficients 

for  A  ,  m  >  0 
m 


colvunn  matrix  of 
expans Ion  In  0 
see  Fqs.  (3-9) 


coefficients  of  cosine 

^^rm  ^r'm^  ' 


A 

m 


chordwlse  variation  of  sine  component  of 
mNth  harmonic  of  y  ;  see  Eqs.  (1.3) 


Ap^  chordwlse  variation  of  radial  velocity 

Induced  on  shroud  by  sine  component  of 
mNth  harmonic  of  bound  blade  vortices; 
see  Eqs.  (l.ll) 


Ap/  chordwlse  variation  of  radial  velocity 

Induced  on  shroud  by  sine  component  of 
mNth  harmonic  of  shed  blade  vortices; 
see  Eqs.  (l.ll) 


Glauert  coefficients  for  B  ,  m  >  C  ; 
see  Eqs.  (3*3) 


(b  )  column  matrix  of  Glauert  coefficients 

B  ,  m  >  0 
m 


(bp^]  column  matrix  of  coefficients  of  cosine 

expansion  In  <t>  of  Bp'  ,  m  >  0  ; 
see  Eqs.  (3.9) 


B 

m 


chordwlse  variation  of  cosine  component 
of  mNth  harmonic  of  y  ;  see  Eqs.  (1.3) 


V 


chordwise  variation  of  radial  velocity 
induced  on  shroud  by  cosine  component 
of  mNth  harmonic  of  shed  blade  vortices; 
see  Eqs.  (1.11) 


c 


shroud  chord;  see  Fig.  1.1 


c 


00 


(j) 


jth  iterate  of  Glauert  coefficients  for 
Cq,  j=0,  1,  2,  see  Eq.  (2.5) 


c 

P 


] 


net  average  pressure  on  shroud,  positive 
radially  outward;  see  Eq.  (2.35) 


colvuun  matrix  of  Glauert  coefficients 
of 


column  matrix  of  jth  iterate  of  Glauert 
coefficients  for  C«,  j=0,  1,  2... 
see  Eq.  (2.20)  ^ 


C 

m 


chordwise  variation  of  mNth  complex 
Fourier  harmonic  of  7  ;  see  Eq.  (1.2) 


C 


0 


(j) 


jth  iterate  of  ;  see  Eq.  (2.4) 


g 


niN 


^mN 


Fourier  integrals  arising  from  shed 
vorticity;  see  Eqs.  (1.12) 


i 


-1  and  dummy  summation  index 


unit  vector  in  r-direction;  see  Fig.  1.1 

"particle  displacements"  for  higher  har¬ 
monic  equations  in  decoupled  form 


VI 


integral  matrix  operator  for  singular 
part  of  ,  identified  as  infinite 
unit  matrix;  see  Eqs.  (3.3) 


iteration  number  and  dummy  svimmation 
index 


propeller  advance  ratio,  u/f2Rp 

integral  matrix  operator  for  singular 
part  of  ;  see  Eqs.  (3*7)  and  (3.8) 

dummy  summation  index 


kernel  function  of  final  governing 
equations  -  for  m  =  C  only  K- 
appears;  see  Eqs.  (1.7)  and  (l  14) 


regular  part  of  K  ;  see  Eqs.  (1.9) 
and  (3.4)  ^ 


regular  part  of  ,  equal  to 


dummy  summation  index 


harmonic  solution  index,  generally 

m  =  0,  1,  2  ... 


matrix  of  eigenvectors  of  [p]  :  see 
Eq .  f  2 . 24  ) 


number  of  blades 


Fourier  coefficients  for  expansion  of 

for  zeroth  harmonic  ~  /^  a  - - 

Eq.  (2.7) 


(^) 


see 


propeller  "forcing  functions"  for  higher 
harmonic  equations  in  decoupled  form; 
see  Eqs.  (3*l) 


elements  of  [p]  ,  identified  as  "curva¬ 
ture  coefficients"  for  m  =  0  ;  see 
Eqs.  (2.10)  and  (2.18) 


integral  matrix  operator  for  regular  part 
of  ;  see  Eqs.  (2.21)  and  (3-5) 

total  induced  velocity;  see  Eq.  (1.1) 


Legendre  function  of  second  kind  emd 
half  order 


integral  matrix  operator  for  regular  part 
of  ;  see  Eqs.  (3*7) 


r/R 

P 

mean  shroud  radius  at  propeller  plane; 
see  Fig.  1.1 


propeller  tip  radius;  see  Fig.  1.1 


^*^mN+i}  '^mN-3/2^ 


forward  flight  speed;  see  Fig.  1.1 


cylindrical  coordinates  fixed  in  pro¬ 
peller;  see  Fig.  1.1 


IC-. 


X  position  of  propeller  pleme  measured 

^  downstream  from  shroud  midchord;  see 

Pig.  1.1 


X  ,r  ,6 
s  s'  s 


cylindrical  coordinates  for  point  on 
shroud  surface 


cylindrical  coordinates  for  position 
of  vortex  element 


X 


x/R  ,  . . . 


a 


n,K' 


coefficients  of  small  argument  expansion 
of  ?  see  Eq.  (1.8) 


y 


strength  of  bound  shroud  vortices; 
see  Fig.  1.1 


r 


strength  of  bound  blade  vortex;  see 
Fig.  1.1 


r/FpU 
("s  - 


C 


(x  -  X  ) 

'  s  p' 

■  ^P  ■ 

local  chordwise  shroud  slope  or  camber; 
see  Fig.  1.1 


€  effective  shroud  camber 

e 


IX 


dununy  sunnnation  index 


shroud  chord  to  diameter  ratio,  c/2B 


eigenvalues  of  [p]  ,  Ix  1  >1x1  >  ... 
see  Eq.  (2.29)  ^ 


successive  approximations  to  X 
see  Eqs.  (2.30) 


diagonal  matrix  with  elements 
Xq  .  Xi  ,  X^  . . .  ;  see  Eq.  (2.25) 

ratio  of  propeller  tip  to  mean  shroud 
radius,  Rp/R 

summation  index  for  Glauert  series, 

V  =  1,  2.  3  ... 


argximent  of  Legendre  function 

(1  +  Ax^V2) 

(1  +  [AXp'*^  +  (1-r^)^^  / 

(1  +  +  (1-?^)^]  / 


"natural  frequency"  for  higher  harmonic 
equations  in  decoupled  form.  mN/J4 


dummy  time  varicd>le 


Glauert  variable,  x  »  -X  cos  <t>  ; 
see  Fig.  2.1 


column  matrix  of  caiia>ar  alsnonta  for  two- 
dimensional  airfoil  with  unit  Glauert 
coefficients;  see  Eq.  (3.6) 


2uigular  rotational  speed  of  propeller; 
see  Fig.  1.1 


denotes  total  differentiation  with 
respect  to  indicated  argument 

transpose  of  colvunn  matrix  or  row  matrix 


Inverse  square  matrix 


GENERAL  HARMONIC  SOLUTIONS 
FOR  THE 

ducted  propeller 

INTRODUCTION 

In  an  effort  to  achieve  a  general  understanding  of 
the  ducted  propeller,  we  have  previously  formulated  a 
three-dimensional  theory^  for  the  forward  flight  regime  at 
zero  angle  of  attack.  The  appropriate  mathematical  model 
was  based  on  the  classical  concepts  of  vortex  propeller 
theory  for  finite  blading  and  of  lifting  surface  theory 
for  thin  wings  in  inviscid,  Incompressible  flow. 

Primary  emphasis  was  placed  on  the  inverse-direct 
problem,  that  is,  the  determination  of  the  shroud  loading 
for  a  given  blade  circulation  distribution.  By  Fourier 
analysis  of  the  radial  component  of  the  propeller  velocity 
field  and  the  bound  vortex  distribution  of  the  shroud,  we 
reduced  the  problem  for  each  harmonic  to  the  inversion  of 
a  coupled  pair  of  line  Integrals  of  the  Cauchy  type. 

The  zeroth  harmonic  corresponds  to  the  representation 
of  the  propeller  as  a  generalized  actuator  disk  having  the 
same  radial  disk  loading  as  its  finite  bladed  counterpart. 
As  such,  it  determines  the  average  shroud  load  or  the  most 
significant  contribution.  The  higher  harmonics  give  the 
time -dependent  loading  on  the  shroud. 
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For  the  zero  harmonic,  the  pair  of  integral  equations 
degenerates  to  a  single  Integral  equation  identical  to  the 


2 


equation  for  an  equivalent  ring  wing  with  a  different 
effective  camber  in  a  uniform  flow.  The  higher  harmonics 
may  also  be  decoupled  into  a  similar  form,  but  with  a 

different  kernel  because  of  the  swirl  of  the  shroud  shed 

2 

vortices.  Since  J.  Weissinger  had  formulated  the  exact 
solutions  to  both  of  these  equations,  the  solutions  were 
left  at  that  time  implicitly  in  terms  of  his  results. 

The  objective  of  the  present  investigation  is  to  ob¬ 
tain  explicit  solutions.  In  particular,  it  is  desirable 
to  find  the  effects  of  the  blade  number  and  loading,  the 
advance  ratio,  the  position  of  the  propeller  plane,  the 
tip  clearance  and  the  shroud  chord  to  diameter  ratio  in 
reasonably  simple  terms.  Except  for  the  influence  of  the 
chord  to  diameter  ratio,  this  is  possible  for  the  zeroth 
harmonic  by  a  study  of  the  effective  Ccunber.  As  the  chord 
to  dicimeter  ratio  appears  in  the  limits  of  the  integral 
in  the  integral  equation,  its  effect  requires  the  inversion 
of  the  equation.  Unlike  the  two-dimensional  thin  air¬ 
foil  problem,  the  Glauert  coefficients  of  the  direct  inver¬ 
sion  of  the  ring  wing  equation  are  all  coupled  by  an  infinite 

2 

set  of  linear  equations  of  infinite  extent  .  Therefore,  we 
attempted  to  achieve  a  simplification  through  cin  iteration 
of  the  corresponding  two-dimensional,  sectional  coefficients 
which  are  independently  determined.  The  saime  would  be  true 
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for  the  higher  harmonics  in  light  of*  the  similarity  of  the 

n 

decoupled  form . 

e 

Here  we  have  set  forth  the  iteration  procedure  in  de¬ 
tail,  as  well  as  an  alternative  method  for  the  higher  har¬ 
monics.  Chapter  One  outlines  a  brief  review  of  the  de¬ 
velopment  of  the  underlying  theory^  and  presents  a  general 
discussion  of  the  basic  equations  to  be  solved.  In  Chapter 
Two,  the  solution  by  iteration  is  found  for  the  zeroth  har¬ 
monic.  Any  desired  order  of  iteration  is  carried  out  in 
an  explicit  form,  and  in  turn,  is  expressed  as  the  expansion 
of  a  single  matrix  operator  of  infinite  order.  This  permits 
the  evaluation  of  the  condition  and  speed  of  convergence 
by  meams  of  the  maotimum  eigenvalue  of  the  matrix  and  the 
identification  of  the  infinite  iterate  with  Weissinger's 
solution.  A  numerical  exeunple  is  given,  the  calculations 
being  based  upon  the  appropriate  included  taibles  of  cur¬ 
vature  coefficients.  These  coefficients  are  related  to 
the  coefficients  tabulated  by  the  Bureau  Technique  Zborowski-", 
but  to  check,  they  are  computed  independently.  Finally, 
Chapter  Three  contains  the  iteration  solution  for  the  higher 
harmonics,  a  direct  solution  by  simultaneous  matrix  in¬ 
version,  and  a  comparison  of  the  features  of  these  solutions. 


CHAPTER  ONE 


BASIC  FORMULATION 


1.1  Mathematical  Model  and  Background 

We  consider  a  ducted  propeller  in  steady  forward  flight 
U  at  constant  rotational  speed  Q  and  zero  incidence  in  an 
inviscid,  incompressible  fluid  otherwise  at  rest.  The  shroud 
and  propeller  blades  are  assumed  to  be  infinitely  thin  and 
the  hub  radius  infinitely  small.  These  assumptions,  together 
with  the  restrictions  of  small  shroud  camber,  moderate  chord 
to  diameter  ratio  and  low  disk  loading  permit  the  conven¬ 
tional  approximations  of  linearized  perturbations . 

In  a  propeller-fixed  coordinate  syst«_ia  x,r,e  ,  the  in¬ 
coming  stream  has  a  uniform  translation  U  in  the  positive 
x-direction  and  a  rotational  component  fir  in  the  positive 
6 -direction,  see  Fig.  1.1.  The  blades  are  at  rest  in  this 
system,  but  the  duct  is  rotating  about  the  propeller  axis 
with  the  angular  speed  fi  .  However,  we  may  set  the  duct 
rotation  equal  to  zero  without  altering  the  potential  solu¬ 
tion  since  we  have  limited  ourselves  to  axisymmetric  ducts 
at  zero  eungle  of  attack^,  consequently,  the  determination 
of  the  flow  is  reduced  to  a  steady  problem. 

The  einalysis  is  formulated  by  the  use  of  appropriate 
vortex  distributions  which  satisfy  the  physical  boundary 
conditions.  For  the  propeller,  we  have  the  classical 
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representation  of  a  bound  radial  vortex  line  of  strength  r(r^) 
for  each  blade  and  the  associated  helical  sheet  of  strength 
-dr/dr^  shed  by  the  line  Where  r^  is  the  radial  distance 
to  any  line  element.  Consistent  with  our  assumptions,  the 
helical  path  of  the  shed  vortex  is  determined  by  the  incoming 
stream.  On  the  other  hand,  the  shroud  is  composed  of  bound 
ring  vortices  7  which  are  concentric  to  the  x-axis  and 
have  a  radius  equal  to  the  mean  duct  radius  R  .  Since  the 
disturbance  field  introduced  by  the  propeller  is  not  uniform 
in  the  angular  direction,  these  bound  vortices  are  accompanied 
by  shed  vortex  cylinders  that  are  composed  of  helical  vortex 
elements  of  strength  .  Analogous  to  the  vortices 

shed  from  the  propeller,  their  helical  pitch  is  fixed  by  the 
incoming  stream. 

With  these  vortex  distributions,  we  may  determine  the 
basic  equation  for  the  shroud  vorticity  in  terms  of  the  pro¬ 
peller  circulation.  That  is,  we  require  that  the  flow  be 
tangent  everywhere  on  the  shroud  surface  and  satisfy  the 
Kutta-Joukowski  condition  at  the  trailing  edge.  With  our 
assumptions,  this  reduces  to  a  feuniliar  thin  airfoil  approxi¬ 
mation, 

U€(Xg)  =  q^  (Xg.R.e^)  .  i^(0g)  (1.1) 

where  €  is  the  local  inclination  angle  of  the  shroud  sur¬ 
face  to  the  x-axis ;  q^  is  the  sum  of  the  velocities  induced 
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by  the  four  vortex  distributions  on  the  mean  shroud  surface 
1  local  unit  vector  in  the  radial 

direction.  The  total  induced  velocity  may  be  determined  by 
means  of  the  Biot-Savart  law  and  Eq.  (l.l)  reduced  to  the 
form  of  an  integral  equation  for  7  with  the  propeller  con¬ 
tributions  transposed  to  the  camber  term. 

1.2  Reduction  of  the  Governing  Equation 

Generally,  the  inversion  of  a  surface  integral  equation 
is  inherently  complex.  However,  in  the  present  case  as  for 
the  ring  wing,  the  closing  of  the  shroud  surface  introduces 
an  essential  simplification.  In  particular,  7(x,0)  may 
be  expanded  in  a  complex  Fourier  series,  or 

00 

7(x,9)su  ^  C„(x)  (1.2) 

m=-« 

where  vanishes  at  the  shroud  trailing  edge  x  =  c/2  and 

N  is  the  number  of  propeller  blades  of  radius  R^  .  Only 
the  harmonic  components  mN  are  necessary  because  of  the 
periodicity  27r/N  of  the  propeller  field.  Since  7  is 
real,  the  complex  coefficients  for  positive  (+m)  and  nega¬ 
tive  (-m)  integers  are  complex  conjugates,  or 


2C  s  B  +  iA 
-m  m  m 


2C  s  B  -  iA 
+m  m  m 


(1.3) 
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UA^  and  UB^  are  the  real  coefficients  of  the  sine  and  cosine 
m  m 

terns  respectively  of  a  trigonometric  Fourier  series  for  7 
over  positive  m  . 

For  m  =  0  ,  we  note  =  0  and  2Cq  =  Bq  .  This 
is  the  zeroth  Fourier  component  of  the  shroud  vorticity  in 
the  propeller-fixed  coordinates,  or  equivalently,  the  time- 
independent  part  in  the  shroud-fixed  coordinates.  As  such, 
its  physical  role  has  been  shown  to  determine  the  time  aver¬ 
age  difference  in  shroud  inner  and  outer  static  pressures. 

Substitution  of  Eq.  (1.2)  into  Eq.  (l.l)  with  the  ap¬ 
propriate  velocity  contributions  permits  the  explicit  integratioi 
over  the  angular  coordinate  and  reduces  the  shroud  contribu¬ 
tions  to  their  Fourier  components  in  .  At  the  same 

time  the  propeller  and  shroud  camber  contributions  may  be 
similarly  reduced.  As  a  result,  the  coefficients  of  these 
Fourier  components  may  be  equated,  forming  in  general,  a  pair 
of  line  integral  equations  for  A  and  B  .  In  other  words, 
our  original  surface  integral  equation  is  replaced  by  a  semi¬ 
infinite  set  of  coupled,  line  integral  equations  for  each 
harmonic^ . 

1 . 3  The  Zeroth  Harmonic  Equation 

Since  A^  =  C  ,  the  integral  equation  for  the  zeroth  har- 
monic  is  particularly  simple.  With  x  =  x  /R  ...  ,  we  have 

o  S 

X 

€(x)=-  /  C-(x)K'(Ax)dx 

e's'  J  C'v'  C'  v'  v 

-X 


(i.M 
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where  €  is  the  "effective"  local  shroud  slope  or  camber; 
k'  ,  the  kernel  for  the  zeroth  haunnonic;  ,  the  sepeura- 
tion  of  the  vortex  element  and  field  point  Ax^  =  (x^-x^)  ; 
and  X  ,  the  shroud  chord  to  diameter  ratio  or  X  =  c/2R  . 

The  effective  camber  is  composed  of  the  physical  shroud 
camber  and  the  axisyrametric  part  of  the  radial  field  induced 
by  the  propeller  at  the  shroud,  or 


€ 

e 


=  € 


(1.5) 


The  propeller  term, 

^  °  sA  -Jq  '  v'  >s'  3'  V  V 

Sj  2  1  +  [AJ2  +  (l-r^)2)/2J^  (1.6) 


depends  on  the  blade  number,  the  advance  ratio  J  =  U/fiR  , 

P 

the  propeller  tip  to  shroud  radius  ratio  [i  ^  R  /R  ,  the 

P 

slope  of  the  blade  circulation  distribution  r'=  d(r/RpU)/dr^  , 
and  the  propeller  plane  position  x^  .  For  convenience,  F 
can  be  expanded  in  a  suitadile  sine  series  and  ex¬ 

pressed  in  terms  of  certain  characteristic  functions.  These 
functions,  involving  numerical  integration^  over  the  Legendre 
function  of  second  kind  and  first  half  order  ,  have  been 

tabulated  over  the  reuige  of  values  of  C.u^x  ^l.C  and 

P 

0.75  <  4  <  0.99  . 
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For  m  =  0  ,  the  kernel  is  simplified  by  the  absence 
of  terms  corresponding  to  helical  vortices  shed  from  the 
shroud.  In  fact,  -k'  is  the  radial  velocity  distribution 
induced  by  a  ring  vortex  of  positive  unit  strength  and  unit 
radius  on  a  cylinder  of  which  the  ring  is  an  element.  In 
terms  of  Qj^  ,  we  have 

■'o  =■  A 


S  1  +  /2  (1.7) 

where  as  in  T'  ,  the  prime  denotes  differentiation  with 

respect  to  the  argxjment  of  the  functions .  as  we  approach 

-  4  S 

the  ring,  dx^  -•  C  and  1  .  But  near  unit  argument  , 

00  00 

Vl,'”)  '  I  “n.K  ®n,K 

K=0  K»0 

For  Eq.  (1.7)  with  n  =  1  we  have  q  =  (5in2  -  4)/2  , 

®1,1  "  (15in2  -  2)/16  ,  ...  and  =  -1/2  ,  3^^  =  -3/16  , 

...  .  Therefore,  the  singular  part  of  -K^  is  equal  to 
l/27rdx^  ,  i.e.,  in  the  immediate  vicinity  of  the  vortex,  the 
flow  field  is  equivalent  to  a  two-dimensional  line  vortex 
tangent  to  the  ring  and  the  integral  of  Eq.  (1.4)  must  be 
interpreted  as  the  Cauchy  principal  value.  On  the  other 
hand,  the  regular  part  of  the  kernel,  say  k'  ,  where 
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-  Kq  +  l/27rAx^ 


(1.9) 


represents  the  three-dimensional  effect  of  the  ring  which 

vanishes  as  the  ring  radius  become  infinite.  behaves 

as  Ax  ^n  Ax  near  the  vortex,  increases  to  a  maximum  ana 

V  V 

then  decays  monotonically  to  zero  as  Ax^  -♦  oo  .  By  inspection 
of  Eqs.  (1.7)  and  (1.8),  both  the  singular  and  regular  parts 
of  the  kernel  are  antisymmetric  about  Ax^  =  C  .  If  we 
compare  Eq.  (1.4)  with  Weissinger's  ring  wing  equation  for 
zero  angle  of  attack,  we  find^  that  the  forms  are  identical 

with  "U^(t])"  equal  to  -27tR'  . 

U  U 

1.4  The  Higher  Harmonic  Equations 

For  m  ,  J  /  0  ,  the  integral  equations  for  and 

are  considerably  more  complex.  The  essential  complica¬ 
tions  arise  from  the  additional  terms  introduced  by  the 
helical  shed  vortices  which  couple  the  integral  equacions,  or 


A.  A 

Ap  +Ap'  =  /  A  K'dx  +0 

Pm  Pm  Jmmv  J 

-X  -X 


B  K  dx 

mm  V 


B, 


A  A 

^  =  B  K'^dx  -  a/  A 

m  Jmmv  J  > 

-X  -X 


K  dx 

mm  V 


(1.10) 


with  a  H  mN/JM.  .  The  term  UAp^  is  the  coefficient  of  the 
mN^^  harmonic  in  the  Fourier  sine  expansion  of  the  radial 
velocity  component  induced  cn  the  duct  by  the  bound  radial 
vortex  lines  which  represent  the  propeller, 


and 
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4*1^ 

UBp '  are  the  coefficients  of  the  mN  sine  and  cosine 
harmonics,  respectively,  in  the  Fourier  expansion  of  the  ra¬ 
dial  velocity  component  induced  on  the  duct  by  the  helical 
vortices  shed  from  the  blades .  These  terms  are  given  by 


Tm 


'm 


I  ri 


Ax  /  — ^  Q  1.  dr 

07,2  P  J  z  3/?  mn-h  V 

^  0  V 


B 


r  m 


I-  f  r 

r2T  J  V 

zh 


— I  r"^  ( S  -  g  )  dr 

2^^2j  j  V  '  mN  ^mN'  v 


(1.11) 


where  the  argument  of  is  of  Eq.  (1.6), 

^mN-^s  general  Legendre  function  of  second  kind  and 

half  order,  and  is  the  sum  of  and  • 

The  functions  g  „  and  h  „  are  Fourier  integrals  defined  by 

mN  mN 


g  „(Ax  ,r  ) 

^mN'  p  v' 

h  „ ( Ax  , r  ) 
mN'  p  v' 

(sin  mNT) 
(cos  mNT) 


mN 


w 

f  [ 


n-r2  2 

2J  u 


®mN-H  ^ 


mN 


(1.12) 


Umjj.^j  and  have  the  argument  cu^  which  is  equal  to  cd^ 

with  AXp  replaced  by  the  quantity  ^  (^^p  “  Ut/r)  ,  the 

Q  „  ,  term  arising  from  the  axial  elements  of  the  helical 
n\N“*5 
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vortices  and  the  term,  from  the  tangential,  t  Is  a 

variable  representing  the  time  for  a  vortex  element,  con~ 

vected  by  the  free  stream,  to  travel  from  the  point  or 
shedding  to  its  position  on  the  helix  and  t  =  ftr  .  For 
|i  ,  and  so  r^  ,  less  than  unity,  there  is  no  question  about 
the  behavior  of  g^^  and  .  When  p.  is  equal  to  unity, 

Eq.  (1.8)  assures  that  the  singular  integrand  is  integrable, 

c; 

amd  the  asymptotic  form"^  of  the  Legendre  function,  or 


TT  (2n):  _ 1_  . 

(n;)2 


(1.13) 


that  it  is  convergent  far  downstream. 

The  kernel  function  may  also  be  defined  in  terms  of 


K  (Ax  )  = 
m'  v' 


_L 

Utt 


g  „( Ax  , 1 ) ] 
^mN'  v  ' 


(1.14) 


where  g  „  has  been  evaluated  at  r  =1  or  on  the  shroud. 
niN  v 

To  understand  the  form  of  Eqs.  (I.IC),  further  descrip¬ 
tion  of  the  kernel  function  is  necessary.  The  term 

-Ax  S  ^.,/47r  is  identified  as  the  mN^^  harmonic  of  the  radial 
velocity  induced  on  the  shroud  by  the  shroud  bound  vortices 
in  non-dimensional  form.  In  other  words,  the  radial  vel¬ 
ocity  of  a  ring  vortex  of  strength  UA^  sin  mN0^dx^  is 
(UA  sin  mN0  dx  )  (-Ax  S*„/47t)  ,  or  "in-phase"  with  the 


strength  variation. 
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In  the  neighborhood  of  the  ring,  the  vortex  again  appears 
as  a  straight,  infinite  line  vortex  tangent  to  the  ring 
with  a  strength  equal  to  the  local  strength,  see  Eq.  (1.8), 
and  the  integrals  of  Eqs.  (I.IO)  are  interpreted  in  the 
Cauchy  principal  value  sense.  As  the  axial  separation  is 
increased,  the  induced  radial  velocity  decays,  the  asymp¬ 
totic  behavior  determined  by  Eq.  (1.13). 

The  other  term  in  and  the  terms  of  K  constitute 

m  m 

the  influence  of  the  shroud  shed  vortices .  Examination  of 
these  terms  reveals  that  the  term  S^^47T  contains  the  "end 
effect"  of  the  semi-infinite  cylinder  of  helical  vortices 
shed  from  our  bound  ring  vortex  at  Ax^  =  0  .  In  other  words, 
as  we  approach  the  ring,  its  shed  vortices  form  an  effective 
planar  sheet  of  vortices  which  have  a  constant  strength  equal 
to  the  rate  of  change  of  the  local  bound  vortex  strength  and 
inclined  at  the  angle  having  the  tangent  1/JM-  to  the  x-axis . 
This  sheet  is  equivalent  to  the  superposition  of  two  semi- 
infinite  sheets  of  vortices,  parallel  eind  perpendicular  to 
the  x-2ocis.  But  by  symmetry  the  parallel  vortices  cannot 
contribute  to  the  "normal"  or  radial  velocity  since  the  con¬ 
tributions  from  either  side  are  equal  and  opposite,  see 
Fig.  1.2.  On  the  other  hand,  the  vortices  perpendicular  to 
the  x-axis  form  an  unbalanced  sheet  and  integrate  to  a 
logarithmic  singularity  given  by  S jjjjj/4'n’  ,  see  Eq.  (1.8). 

The  remaining  terms 

volved  in  Eqs.  (I.IC)  are  both  regular,  the  behavior  of 


c- 
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FIGURE  1.2 


DECOMPOSITION  OF  THE  LOCAL  SHED  VORTEX  SHEET 


following  from  our  previous  discussion  for  the 

propeller  for  zero  tip  clearance.  It  may  be  shown  that 
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mN 


(1.15) 


and  consequently,  from  Eqs.  (1.8),  (1.12)  and  (I.13), 

is  regular.  The  logarithmic  singularities  of 
the  integreuids  of  both  g^(Ax^,l)  and  hj^(Ax^,l)  may  be 
interpreted  analogous  to  S .  That  is,  for  a  field 
point  downstream  of  a  vortex  ring  at  x^  ,  the  shed  vortices 
from  the  ring  locally  appear  as  an  infinite  planar  surface 
of  constant  strength  vortices  inclined  to  the  x-axis.  Only 
this  time,  by  symmetry,  both  the  contributions  to  the  "normal" 
or  radial  velocity  from  the  component  sheets  parallel  and 
perpendicular  to  the  x-axis  cancel . 

The  kernel  function  or,  more  important  as  we  will 

see,  its  derivative  with  g^  given  by  Eq.  (I.15)  has 

not  been  evaluated  in  general  in  closed  form.  For  the  special 
case,  however,  of  infinite  advance  ratio,  the  integral  equa¬ 
tions  decouple,  the  terms  in  K  vanishing.  The  terms 

m  m 

remain  and  h^  reduces  to  simply  an  integral  over  < 

or  a  cylinder  of  semi-infinite  straight  line  vortices  parallel 
to  the  x-axis.  This  may  be  integrated  explicitly^  cind  yields 
the  result, 

limv's-  —  -  1  _  mN 

j-*ai  %  ~  27T  2ttAx^  it 


(1.16) 


Where  is  the  function  introduced  by  J.  Weissinger  in 

his  ring  wing  theory.  It  is  a  regular  function,  increasing 
with  a  logarithmic  slope  from  zero  at  the  origin  to  an 
asymptotic  value  of  7nnN/2  as  -►  »  for  m  >  0  .  As 

m  -►  «  ,  approaches  the  general  form  of  a  unit  step 


function 


CHAPTER  TWO 


DETERMINATION  OF  THE  NET  AVERAGE  DUCT  LOAD 

2.1  General  Considerations 

If  we  neglect  the  regular  part  of  the  kernel  for  the 
zeroth  harmonic,  Eq.  (1.4)  reduces  by  means  of  Eg.  (1.9) 
to  the  integral  equation  of  the  classical  two-dimensional 
thj.n  airfoil  problem.  To  solve,  we  introduce  a  new  axial 
coordinate  <t>  such  that 

X  =  -  X  cos  0  (2.1) 

where  0  <  0  <  tt  ,  see  Fig.  2.1.  Expressing  the  vortex 
distribution  Cq(x)  in  a  Glauert  series,  we  can  then  in¬ 
vert  the  integral  equation  and  obtain  a  solution  explicitly 
in  terras  of  the  Fourier  cosine  coefficients  of  the  effective 
camber  distribution. 

As  we  have  pointed  out,  k'  represents  the  three- 
dimensional  effect  of  the  shroud.  Thus,  we  can  consider 
this  contribution  as  a  correction  to  be  combined  with  the 
effective  camber,  permitting  an  iteration  solution  of 
Eq.  (1.4)  in  the  form  of  the  Glauert  series.  That  is, 
we  may  write  from  Eq.  (2.1) 
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€  +  X 
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/  =0  *0  K 


71 

1-  r 

27r  J 


K 

cos  0  -  cos  0_ 

V  S 


(2.2) 


where  the  regular  part  of  the  kernel  Is 


(2.3) 


and  is  implicitly  a  function  of  0^  .  Then,  if  we  set 

Cq  equal  to  zero  on  the  left  hand  side  of  Eq.  (2.2),  we  can 
find  a  "zeroth"  approximation  from  the  two-dimensional  theory. 
Inserting  this  distribution  back  into  the  left  hand  side,  we 
determine  a  first  approximation  which  includes  a  correction 
due  to  the  surface  curvature.  Subsequently,  the  procedure 
may  be  repeated.  In  the  limit,  it  will  converge  to  the  exact 
solution  of  Eq.  (1.4),  or 


'0 


j-oo  ^0 


(2.4) 


where 


(j  )  =  c  (j  )  I  +  cos  <P 


00 


sin  0 


00 

I 

v=l 


'Ov 


(j) 


sin  v0 


(2.5) 


and  the  superscript  (j)  denotes  the  iteration, 
j  ®  0,  1,  2,  ...  .  Such  a  solution  possesses  the  familiar 
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leading  and  trailing  edge  behavior  of  an  infinitely  thin, 
two-dimensional  airfoil,  l.e.,  a  square-root  singularity  at 
the  leading  edge  and  vanishing  perturbational  velocity  at 
the  trailing  edge. 


2.2  Expansion  of  the  Kernel  Function 

To  facilitate  the  iteration  we  expand  in  the  square 

0  <  0  <  TT  and  0  <  0  <  tt  as 

““  S  —  V  "" 

00  00 

^0  Z 

k=0  i=0 

The  Fourier  coefficients  p^  orthogonality  are 

TT  TT 

P-k  £  \  f  /  ^0  ‘^'^V  ‘^'^S 

0  0 


for  k  ,  e  0  :  and  ^  are  one-half  and 

one-quarter  of  these  values.  Since  is  an  odd  function 

in  Ax^  ,  see  Eqs .  (1.7)  and  (1.9).  we  may  interchange  0^ 
and  0^  in  Eqs.  (2.3)  and  (2.7)  and  find, 


k,£ 


£  ,k 


(2.8) 


Consequently,  p  ,  p.  ,  ,  p^  _  . . .  all  identically  vanish. 
In  addition,  we  may  replace  0^  by  (7T-0g)  and  0^  by 
(7T-0  )  in  Eq.  (2.7)  and  show  that  p,  .  =  0  if  (k+£) 

V  A , 


is  an  even  integer. 
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These  coefficients  may  be  related  to  the  work  of 
2 

J.  Welsslnger  who  expanded  an  Integral  over  a  weighted 

kernel  function  instead  of  the  regular  part  of  the  kernel . 
From  this  integral  designated  by  "fy"  /  we  find  after 
comparison  of  nomenclature. 


P 


k,  i 


(2.9) 


where  the  "b  ’s"  are  the  coefficients  he  used  with  v  =  k 

V  ,  [J. 

and  =  I  for  v,n.  /  0  but  reduced  like  the  p^  ^'s  in 
Eq.  (2.7)  for  vanishing  v  and/or  p.  .  The  Bureau  Technique 
Zborowski  has  tabulated-^  the  "b  for  \  =  0.5»  l.C,  1.5» 

V  ^  p, 

and  2 , C  and  v,p=0,  1,  ...  7  . 

Subsequently,  we  will  find  that  these  coefficients 
as  such  do  not  appear,  but  rather  the  combinations  given  by, 


\,C  ^"Pk,l  ^^k,C^ 


^k,l  .  'Pk,2  ^Pk,C^ 


^k,2  2  ^Pk,3  ^k,l^ 


k,  £+1 


^k, £-1 


) 


(2. 1C) 
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t 


except  for  k  =  0  ,  in  which  case  they  are  of  opposite  sign. 
Their  values  for  k,i  *  0«  1«  •••  7  have  been  computed 
from  Eqs.  (2.7)  and  (2.10)  by  a  releocation  procedure^  carried 
out  xintil  the  difference  between  successive  computations  is 
less  than  +0.00001.  The  results  over  the  range  of  X  =  0.25, 
0.50,  0.75  and  1.0  are  presented  in  Tables  2.1,  2.2,  2.3 
and  2.4  respectively.  For  X  =  0.50  euid  1.0,  they  are  in 
general  the  same  as  found  by  means  of  Eq.  (2.9)  from  the 
tables  of  the  Bureau  Technique  Zborowski  within  a  discrepancy 
of  0(0.001)  or  less,  except  for  a  few  diagonal  terms  at  the 
higher  values  of  k  and  i  . 


2.3  Iterative  Solution 

We  can  now  pursue  the  development  of  the  iterative 
solution  as  previously  outlined. 

First,  we  reduce  Eq.  (2.2)  to  the  two-dimensional  form 
for  the  zeroth  approximation  ,  or 
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sin 

cos  <t>  -  cos  0 

V  s 


d0 
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(2.11) 


7 

The  solution  is  given  by  the  classical  form  of  H.  Glauert  , 
or  Eq.  (2.5)  with  j  =  0  .  To  determine  the  coefficients, 

is  substituted  into  Eq.  (2.11).  After  suitable  mani¬ 
pulation  of  the  integrand  by  trigonometric  identities,  the 
integration  may  be  carried  out  term  by  term  by  means  of  the 


well  known  relation 
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-0.C0003  0  0.00006  0  -0.00330  0  0.00550 
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(2.12) 


for  K=0,  1,  2,  ...  .  As  a  result,  the  contribution  for 
the  singular  component  of  the  vortex  distribution  is  constant; 
i.e.,  the  xiniform  induced  velocity  for  a  flat  plate  in  two- 
dimensional  flow.  The  contribution  for  each  sine  harmonic 
is  a  corresponding  cosine  harmonic.  Or, 
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By  orthogonality,  then. 
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(2.1^) 


for  V  =  1,  2,  ...  .  Eq.  (2.14)  together  with  Eq.  (2.5) 
with  j  =  0  provides  the  zeroth  approximation  with  as 

a  function  of  0  given  by  Eqs .  (1.5),  (1.6)  and  (2.1). 

Note  that  the  signs  of  Eq.  (2.13)  are  reversed  from  the  usual 
two-dimensional  airfoil  result^  because  we  have  defined 
the  positive  sense  of  the  bound  vortices  according  to  the 
right  hand  screw  rule,  see  Fig.  2.1. 
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Next,  we  replace  Cq  by  in  the  left  hwd  side 

of  Bq.  (2.2)  and  obtain  an  equation  for  the  first  approximation 

to  the  boxind  vortex  distribution.  This  equation  is  the 
same  as  Eq.  (2.11)  except  that  the  effective  cauriber  is  changed 
by  the  additional  curvature  term.  To  solve,  we  substitute 
Eq.  (2.5)  with  j  =  1  and  proceed  as  before.  The  right  hand 
side  reduces  correspondingly  to  the  form  of  the  right  hand 
side  of  Eq.  (2.13)  with  the  new  coefficients  Cqq^^^  and 
.  For  the  left  hand  side,  we  use  Eq.  (2.13) 
and  the  expansion  of  Eq.  (2.6)  for  k'  .  The  integration 
over  Kq  is  carried  out  by  means  of  the  integrals. 


t4>  d<t> 

V  V 


fr  , 

0  , 


i  =  0 

I  =  1,  2,  3»  ••• 


(2.15) 


£<p 

V  V 


^/2  , 

0  , 


i  =  1 

i  -  0,  2,  3f 


(2.16) 


TT 

r  , 

j 

0 


V0  COS 
V 


£0  sin  0  d0 

V  V  V 


V2  , 

0  , 


-7r/4  . 


£  =  V  -  1  =  0 
£  =  V  -  1  >  0 
£  V  ±  1 

£=  V  +  1  (2.17) 


and  yields  in  turn  a  cosine  series.  Equating  coefficients 
of  like  terms,  we  have 
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c  =  c  +  P  c  +  P  c  + 

'^OO  ^00  ^0,0  ^00  *^0,1  *^01  ** 


c  -  c  (°)  +  P  c  (0)  +  P  c  + 

01  01  ^1,0  00  ^1,1  ^01  +  ••• 


(2.18) 


Eqs.  (2.18)  reveal  that  the  first  approximation  to  any  co¬ 
efficient  of  our  bound  vortex  distribution  for  the  zeroth 
shroud  harmonic  is  equal  to  the  two-dimensional  value  of 
the  coefficient  plus  a  correction.  The  correction  is  de¬ 
termined  from  all  the  coefficients  of  the  zeroth  approxi¬ 
mation  weighted  by  the  curvature  coefficients  ^  . 

In  other  words,  each  of  the  basic  vortex  distributions 
which  represent  the  effective  camber  in  the  two-dimensional 
case  induce  an  additional  camber  when  introduced  as  a 
three-dimensional  distribution. 

This  process  may  be  repeated  successively.  Since 
only  the  superscript  varies,  we  have  the  expressions. 


c  (j)  =  c  +  P  c  +  P  c  + 

00  00  *^0,0  00  *^C,1  01 


'01 


(j)  =  c  +  P  c  +  P  c  + 


'00 


1,0  oc 


1,1  01 


(2.19) 
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til 

for  the  j  approximation  in  recursion  form.  Eq.  (2.19)  to¬ 
gether  with  Eqs.  (2.4)  and  (2.3)  then  constitute  the  exact 
solution  for  Cq  if  the  limit  exists. 

For  a  symmetric  effective  airfoil  or  €  (-x)  =  €  (x)  ,  we 
have  “*  ^  ®  from  Eqs.  (2.14).  But 

Pq  i»  ^03'  ^05  '**  identically  zero  as  well  as  ^2  1' 

^2  3'  ^2  5  “*  ^4  1'  ^4  3'  ^4  5  so  forth,  see 

Eqs.  (2.10).  From  Eqs.  (2.19)»  then, 

...  =0  for  any  value  of  j  ,  reducing  the  Glauert 
series  of  Eq.  (2.5)  to  only  the  terms  for  v  =  1,  3,  5  •••  • 
Therefore,  the  net  loading  for  the  zeroth  harmonic  would  be 
symmetric  about  the  midchord  and  zero  at  both  the  leading  and 
trailing  edges  as  required  by  consideration  of  the  corresponding 
solution  for  the  equivalent  ring  airfoil  in  reverse  flow. 

2.4  Convergence  of  the  Iteration 

From  physical  considerations  the  convergence  of  the  solution 
depends  on  X  .  If  the  shroud  di2uneter  is  infinitely  large  com¬ 
pared  to  the  shroud  chord,  ^  ^0^^^  "converges" 

to  the  exact,  two-dimensional  solution  in  a  single  step.  As 
X  increases,  the  three-dimensional  effect  of  the  shroud  in¬ 
creases,  reducing  the  domination  of  the  Cauchy  integral  in 
Eq.  (2.2).  The  speed  of  convergence  is  expected  to  diminish 
accordingly  and,  in  fact,  to  vanish  as  X  -►  «>  . 

To  examine  this,  we  use  matrix  operations.  We  define  the 
infinite  column  matrix  of  the  Glauert  coefficients. 
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r  ^ 

(j) 


'00 


'01 


(j) 


'02 


(j) 


(2.20) 


V. 


and  the  infinite  square  matrix  of  the  curvature  coefficients 


^0,0 

^0,1 

p 

0,2 

• 

[P]  H 

"ipO 

^1,1 

^1,2 

^2,0 

^2,1 

^2,2 

I  J  (2.21) 

With  Eqs.  (2.20)  and  (2.21),  Eqs .  (2.19)  become 

If  we  apply  Eq.  (2.22)  to  the  first  and  subsequent  approxima¬ 
tions  and  use  the  unit  matrix  [ of  infinite  order,  we  have 

(Cq<J))  =  (tiJ  +  [Pl  +  [pj®  +  ...  tpjj)  (2.23) 
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That  is,  each  iteration  corresponds  to  a  successive  appli¬ 
cation  of  the  matrix  of  the  curvature  coefficients.  Hence, 
the  iteration  will  converge  if  and  only  if  the  limit  of  the 
sum  of  these  operations  as  j  ,  or  Neumann  series, 

o 

exists  . 

Letting  [m]  be  a  matrix  with  columns  composed  of 
the  eigenvectors  of  [p]  ,  we  have 

[p]  =  [m]  TaJ  [m]"^  (2.24) 


where  [m]”^  denotes  the  inverse  of  [m]  and  TaJ  is 
a  diagonal  matrix. 


rAj 


0 


0  0 


0  0 


(2.25) 


of  the  associated  eigenvalues  of  [p]  or  X^,  X^,  X^  ...  . 
These  elements  are  numbered  such  that  1 X^l  >  1 X^l  >  |  X^l  ...  , 
where  the  eigenvalue  in  [A]  is  in  the  same  column  as  its 
corresponding  eigenvector  in  [P]  . 

By  mecins  of  Eq.  (2.24),  we  may  express  the  repeated 
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operation  of  [p]  on  itself  by, 

[P]^  »  [m]  f-AjJ  [m]"^  (2.26) 

Substitution  of  Eq.  (2.26)  into  the  term  of  Eq.  (2.23) 
then  yields, 

=  [m]  (TiO  +  TaJ  +  +  ... 

...  +  TA-J^)  [M]“^  (2.27) 

with  ( TlJ  +  CAJ  +  ["Aj^  +  ...  rAj-^)  equal  to, 

(1+Xo+Xq+. . .X^ )  0  0 

0  ( 1+X j^+X^+ . .  .X^  )  0 

0  0  (l+X2+x|+. . .x^) 

By  inspection,  the  limit  exists  if  the  limiting  elements 

p  p 

(1+Xq+Xq+.  . . ) ,  ( 1+Xj^+X^+. . . ) ,  ...  exist,  or  equivalently  if 

UqI  <  1 


(2.28) 
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since  is  the  greatest  eigenvalue  in  absolute  magnitude. 

This  inequality  constitutes  a  crlterlm  £or«  and  assures 

the  convergence  of,  the  iteration  process. 

In  principle,  we  can  determine  X^  for  each  chord  to 
diameter  ratio  exactly  from  the  secular  equation  of  [p]  or 


(2.29) 


In  practice,  however,  we  generally  have  to  resort  to  an 
approximate  calculation.  One  method  is  to  assume  a  trial 
eigenvector  of  unit  elements  and  operate  with  [p]  succes¬ 
sively.  If  we  choose  a  corresponding  eigenvalue  X^^^^  of 
unity  and  normalize  the  leading  element  at  each  step,  we  have 
for  the  first  subsequent  approximations. 


X  (^)  =  ry.  D  ^  /  Tk  (0) 
0 


0*"’  =  (f  f  ^0,J 


0 


(1) 


Xo(3)  . 


=  ( 


k  j  i  ^0,k  ^k,j  \ 


(2) 


(2.30) 


We  have  computed  X^  for  X  =  0.25,  0.50,  0.75  and  1.00 
from  Eqs.  (2.30)  and  Tables  2.1,  2.2,  2.3  and  2.4  respectively. 
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The  results  are  given  in  Fig.  2.2  and  Table  2.5  and  overlap 
for  X  «  0.5#  1.0#  1.5#  and  2.0  conqputed  fron  the  BTZ  data^. 

Immediately  we  see  from  Eq.  (2.38)  that  our  solution 
converges,  at  least  for  0  ^  X  ^  2  .  We  see  also  that  the 
speed  of  convergence  decreases  as  anticipated.  That  is, 
from  Eq.  (2.27)  X^-^  determines  the  accuracy  of  the  j 

iteration.  Since  X^  increases  with  increasing  X  ,  the 
nuinber  of  iterations  for  convergence  to  a  prescribed  ac¬ 
curacy  increases . 

We  have  no  numerical  data  beyond  a  chord  to  diameter 
ratio  of  two.  On  the  other  hand,  the  radial  flow  induced 
by  a  vortex  element  of  a  vortex  ring  through  a  point  on  the 
shroud  can  not  be  as  "large"  as  the  "infinite"  contribution 
of  the  element  at  the  point  itself  unless  the  shroud  radius 
is  zero.  Consequently,  the  extremum  eigenvalue  should  con¬ 
tinue  to  increase,  approaching  unity  asymptotically  as 
X  -►  00  .  We  will  see  later  that  X^^  1  .  Physically,  it  can 

not  return  to  zero. 

2.5  Example 

To  illustrate  the  iteration  procedure  we  have  computed 
the  net  shroud  pressure  distribution  of  the  average  load  for 
a  particular  ducted  propeller  with  cin  effective  conical 
Ccunber . 

The  distribution  selected  for  the  propeller  circulation 
is  the  typical  distribution  which  was  previously  used^, 


FOR  ZEROTH  HARMONIC 
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TABLE  2.5 

EXTREMUM  EIGENVALUE  \q  OF  [pj 


X 

1 

A 

i 

"0 

PRESENT 

BTZ 

0.00 

0.000 

0.000 

0.25 

0.030 

0.50 

0.086 

0.086 

0.75 

0.148 

1.00 

0.210 

0.210 

1.50 

0.317 

2.00 

1 

0.402 

.59 


or 


f  -  (2/3>/3)  (2  »in  ni  -  sin  27r£)  (2.31) 

where  r  =  r/R^  .  It  represents  a  disk 
inflow,  of  41.1  psf  at  U  =  200  mph  and 
of  the  parameters  we  take, 

J  =  1 

N  =  2 

=  -  1.0“ 

e 

1  =  2/3 

4  =  0.95  (2.32) 

The  corresponding  geometrical  shroud  camber  €  has  been 
plotted  in  Fig.  2.3-  The  actual  shape  of  the  shroud  can  be 
obtained  by  numerical  integration. 

Only  is  different  from  zero.  From  Eqs .  (2.14) 

and  (2.20),  the  zeroth  approximation  to  the  Glauert  coefficients 


loading,  neglecting 
STP.  For  the  values 


r 
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-0.0349 

0.0000 

0.0000 


(2.33) 


To  find  the  values  of 
X  =  0.50  and  X  =  0.75 


^  ,  we  interpolate  the  data  for 
of  Tables  2.2  and  2.3  and  obtain 


0.104  0 

0.208  0.126 

0.045  0 

0.000  -0.008 

0.000  0 

0.000  0.000 

0.000  0 


0.052  0 

0  -0.023 

0.030  0 

0  0.011 

-0.004  0 

0  -0.002 

0.000  0 


0.000  0 

0  0.000 

-0.008  0 

0  -0.004 

0.006  0 

0  0.004 

-0.001  0 


0.000  .  .  . 

0 

0.000 

0 

-0.002 

0 

0.002 


Now  [Cq)  may  be  approximated  by  means  of  Eq.  (2.23) » 
the  number  of  iterations  determined  consistent  with 
and  the  accuracy  afforded  by  [P]  .  This  requires 


<  0.0005  or  j  =  4  since  =  0.128  ,  see  Fig.  2.2. 
Accordingly,  ([l]  +  [P]  +  ...  [P]^)  i» 


1.119 

0 

O.OoO 

0 

0.000 

0 

0.000 

0.266 

1.144 

0.015 

-0.026 

0 

0.000 

0 

0.052 

0 

1.033 

0 

-0.008 

0 

0.000 

-0.002 

-0.009 

0 

1.011 

0 

-0.004 

0 

0.000 

0 

-0.004 

0 

1.006 

0 

-0.C02 

0.000 

0.000 

0 

-0.002 

0 

1.004 

0 

0.000 

• 

• 

• 

0 

0.000 

0 

-0.001 

0 

1.002 

and 

r 

-0.0391 

-0.0093 

-0.0018 


0.0001 
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These  results  are  substituted  Into  the  Glauert  series 
emd  the  difference  In  the  average  pressure  coefficient  across 
the  shroud  c  ]  ,  taken  positive  radially  outward,  Is  com- 

Ir 

puted  from  the  relation^. 


c 

P 


(2.35) 


It  has  been  plotted  in  Fig.  2.4,  together  with  the  cor¬ 
responding  distribution  for  the  effective  airfoil  in  two- 
dimensional  flow,  i.e.  j  =  0  .  The  influence  of  the  shroud 
curvature  is  about  39%  at  the  midchord. 

In  Table  2.6  we  have  also  presented  c^]  for  each  step 
in  the  iteration  process  to  show  the  speed  of  convergence 
and  to  verify  the  accuracy  of  the  solution. 


2.6  Direct  Inversion 

j.  Weissinger  first  achieved  an  exact  solution  to  the 

P 

ring  airfoil  problem  by  direct  inversion  .  His  solution  may 
be  related  to  our  iteration  procedure  in  the  following  manner . 

When  a  Neumann  series  converges,  we  know  that  the  limit 

0 

is  expressible  as  an  inverse  matrix^,  or 


(Tij  -  [p])’^  =  rij  +  ip]  +  [p]^  +  ... 


(2.36) 


as  j 


00 


,  then,  Eq.  (2.23)  becomes  from  Eqs .  (2.4)  and  (2.36) 
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(oq)  -  (tl)  -  (2.37) 

This  is  the  same  as  Welsslnger's  result  gotten  by  substituting 
a  Glauert  series  simultaneously  into  both  sides  of  Eq.  (2.2) 
and  equating  like  cosine  terms,  cf.  Eq.  (2.19)  without 
the  superscripts  (j)  and  (j-l)  • 

The  form  of  Eq.  (2.37)  permits  two  conclusions; 

(i)  The  difference  in  average  pressure  across  the  shroud  is 
identical  to  that  of  an  equivalent  airfoil  with  a  modified 
camber  in  two-dimensional  flow.  To  derive  the  modified 
camber,  we  determine  the  three-dimensional  coefficients  and 
replace  them  for  the  corresponding  two-dimensional  values 
in  Eq.  (2.13)*  (li)  Since  we  know  the  solution  by  direct 
inversion  exists,  the  inverse  matrix  ([l]  -  [p])""^  must 
exist  and  1  [l]  -  [pJl  0  .  Consequently,  comparison  with 
Eq.  (2.29)  precludes  the  possibility  that  =  1  . 

Either  iteration  or  direct  inversion  may  be  used  in 
amy  application,  each  becoming  progressively  cumbersome 
as  X  increases.  For  numerical  purposes,  iteration  affords 
a  simple  evaluation  of  the  computational  accuracy  amd  minimizes 


the  chances  for  arithmetical  errors . 


CHAPTER  THREE 


HIGHER  HARMONIC  SOLUTIONS 


3.1  Methods  of  Approach 

Neither  the  iteration  process  nor  the  equivalent  direct 
inversion  of  Eq.  (1.4)  for  the  zeroth  harmonic  depends  upon 
the  nature  of  the  regular  part  of  the  kernel.  Therefore, 
if  we  can  decouple  Eqs.  (I.IO)  for  the  higher  harmonics 
into  a  form  with  a  kernel  which  has  a  similar  Cauchy  sin¬ 
gularity,  we  can  utilize  the  same  technique.  Fortunately, 
we  found  such  a  way  to  decouple  these  equations  earlier^ 
and  so  have  a  solution  in  effect  for  the  higher  harmonics. 

An  alternate  method  of  solution  is  possible  by  direct 
inversion  of  the  coupled  integral  equations.  That  is,  we 
follow  the  development  of  the  integral  of  the  product  of 
the  Glauert  series  and  the  kernel  function  itself,  as  well 
as  its  derivative,  into  a  cosine  series  in  0  and  then,  de¬ 
termine  sirault2meously  the  coefficients  of  like  terms. 

Generally  speaking,  either  of  these  approaches  is 
quite  complicated.  One  simplification  is  achieved,  though, 
if  we  take  the  limit  of  infinite  advance  ratio.  In  this  case 
Bp,  B^,  ...  =  C  and  the  integral  equation  for  Aj^, 

Ap,  A^,  ...  is  the  same  as  for  an  equivalent  asymmetric 
ring  wing  at  zero  angle  of  attack  in  uniform  flow,  the 


effective  asymmetry  arising  from  the  circumferential 
variation  of  radial  velocity  induced  by  the  propeller. 
Accordingly,  such  a  limit  is  of  particular  interest  to 
assess  the  importance  of  the  higher  harmonics  and  to  check 
the  computational  procedures  of  numerical  programs . 

3.2  Solution  by  Reduction  to  Zeroth  Harmonic  Form 

We  will  only  outline  here  the  reduction  of  the 
coupled  equations  to  the  zeroth  harmonic  form . 

•  We  let  the  second  integrals  on  the  right  hand  side 

of  Eqs.  (1.10)  be  denoted  by  I  and  I  respectively. 

£\  D 

Since  the  derivative  of  Kj^(Ax^)  with  respect  to  Ax^ 

is  equal  to  the  derivative  with  respect  to  x^  which 

can  be  brought  out  from  under  the  integral  sign  of  the 

first  integrals,  Eqs.  (I.IO)  so  represent  a  set  of  first 

order,  linear  differential  equations  in  x^  .  To  separate 

,  we  take  the  derivative  of  the  first  equation  and 

subtract  o  times  the  second;  or  to  separate  I  ,  we 

B 

multiply  the  first  by  a  and  add  to  the  derivative  of 

the  second.  Both  results  are  of  the  form  of  a  second 

order,  linear  differential  equation  for  a  simple  harmonic 

oscillator  subject  to  a  forcing  function.  That  is,  we 

view  I  and  I  as  "particle  displacements"  and  the 

axial  coordinate  x  as  "time".  The  "motion"  is  un- 

s 

damped  with  a  "natural  frequency"  o  .  The  "forcing 
functions"  are, 
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'’b  ^  m)'  ^  (^rm  ^  Ap'n,)  (3.1) 

From  these  differential  equations,  solutions  for 

and  I  may  be  found  separately  which  provide  decoupled 

equations  for  A„  and  . 

M  M 

For  the  inhomogeneous  solutions .  we  determine  the 

"indicial  response"  to  a  unit  step  function  and  use  Duharoel's 

method  of  superposition.  If  the  resulting  equations  are 

differentiated  with  respect  to  x  and  the  differentiation 

s 

is  carried  back  inside  the  integrals  and  ,  the 

final  equations  assume  the  form, 


X 

r  ® 

X 

r 

/ 

cos  Ox 

dx  = 

I 

A  (x  ) 
m'  V ' 

-X 

-X 

h 

X 

r 

J 

PpU) 

cos  ax 

dx  = 

j 

3  (x  ) 
m'  V ' 

-X 

A 

K^i,Ax  )  dx 
m'  v'  V 


X  )  dx 

V '  V 


(3.2) 


Eqs .  (3.21  are  now  identical  in  structure  to  Eq.  (1.^) 
and  yield  to  solution  as  before,  though  in  detail  they 
are  quite  complex. 

comparing  the  terms  on  the  left  hand  side  by  means  of 
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Eqs.  (1.5)»  (1*6),  (1.11)  emd  (3*1) »  we  see  that  several 
additional  terms  as  well  as  two  Integrations,  i.e.  on  t 
2md  X  ,  are  Involved.  These  integrations  are  further 
aggravated  by  the  "integral"  presence  of  the  advance  ratio 
J  and  the  blade  nviihber  N  .  In  the  zero  harmonic,  they 
appear  simply  as  factors  of  the  propeller  contribution  to 
€ 

e 

On  the  other  hand,  the  relative  complexity  of  the 

terms  on  the  right  hand  side  of  Eqs.  (1.4)  and  (3*2)  is 

not  so  great.  Only  a  few  more  terms  and  an  Integration  on 

T  are  imposed.  Previously  no  parameters  appeared  in  the 

kernel.  Now  we  have  j  ,  p.  ,  m  and  N  or  explicitly,  the 

combinations  JU  and  mN  ,  which  along  with  X  enter 

in  the  curvature  coefficients.  We  also  notice  that  the 

kernel  is  not  antisymmetric  about  Ax^  =  0  because 

of  g  ,1)  or  h  „(Ax  ,1)  ,  see  Eq.  (1.12).  This 

extends  the  calculation  of  the  corresponding  Fourier 

coefficients  p^  ^  for  the  expansion  of  the  regular  part 

of  the  kernel  k'  to  include  all  values  of  k  and  I  , 

m 

no  coefficients  vanishing  identically. 

Besides  the  inhomogeneous  solutions  to  the  differential 

equations  for  and  ,  there  are  also  homogeneous 

solutions.  The  corresponding  load  distributions  which  they 

produce  may  be  found  in  the  same  way  as  the  A  amd 

mm 

for  the  inhomogeneous  solutions,  but  the  physical  interpre¬ 
tation  of  these  distributions  is  not  yet  clear. 
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3.3  Direct.  Simultaneous  Inversion 

Rather  than  decouple  the  governing  equations  by  differ > 

entiation,  we  ceui  express  them  in  matrix  form  as  for  the 

zero  harmonic  and  solve  by  direct,  simultaneous  inversion,  i.e. 

we  assume  both  A  and  B  in  the  form  of  a  Glauert  series 

m  m 


A  =  a 


m  mO  sin  <t> 


00 

i  ^  +  y  a  „  Sin 

Sin  <p  Z_i  mV 


v=l 


B  = 
m 


00 

1  ..cos  »  ^  y  b  „  Sin 
mO  sin  0  Z_i  rov 
v=l 


(3.3) 


and  substitute  in  Eqs.  (1.10).  With  suitable  Fourier  ex¬ 
pansions  of  K  and  k'  the  integrals  over  A„  and 

m  ra  ^  mm 

may  be  written  as  pure  cosine  series  in  the  Glauert  angular 

variable  0  for  the  position  on  the  shroud.  The  propeller 
s 

terms  An  ,  An'  and  Bn'  can  be  expressed  in  a  similar 
1  m  1  m  1  m 

fashion.  Equating  like  coefficients,  we  may  write  the  set 
of  linear  equations  in  the  Glauert  coefficients  as  a  matrix 
equation  and  invert. 

As  noted  earlier,  K  has  a  logarithmic  singularity, 

m 

see  of  Eq.  (1.14)  and  Eq.  (1.8).  Consequently,  if  we 

define  a  regular  part  , 


2 


(3.^0 


5? 


the  terms  In  k'  may  be  reduced  right  away  from  the 

XR 

zeroth  harmonic,  or  from  Eqs.  (2.1),  (2.12),  (2.15), 

(2.16),  (2.17)  together  with  Eqs.  (3.3)  and  R'  ex- 

m 

pressed  in  the  double  Fourier  series  of  Eq.  (2.6),  cf. 

Eqs.  (2.13)  and  (2.37), 

X 

Ta  k'  dx  =  -  i(j> [o^-]  {a  ) 
j  ra  m  V  s  K  m 

-X 

X 

/  B  K  dx  =  -  {0  )'^  [0„<]  [b  } 

J  m  m  V  s  K 

-X 

[Oj^/]  =  [“lO  -  [p]  (3.5) 


T 

The  matrix  (0  }  is  defined  by, 

9 


r 


-cos  0 


^‘^s^  2  1  -cos  2  0. 


(3.6) 
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where  "T"  denotes  a  transpose  or  row  matrix.  is 

the  \anit  matrix  and  [p]  is  the  matrix  of  Eq.  (2.21)  having 

the  elements  of  Eqs.  (2.?)  and  (2.10)  with  in  place 

of  .  The  remaining  matrices  (a  }  and  (b  }  are  column 
u  mm 

matrices  of  the  form  of  Eq.  (2.20)  with  the  Glauert 
coefficients  of  Eqs.  (3*3)  respectively  as  elements. 

For  the  reduction  of  the  terms  in  of  Eqs.  (I.IO), 

we  proceed  in  a  similar  manner  and  separate  the  singular  and 
regular  parts.  The  logarithmic  term  arises  in  the  calculation 
of  the  average  shroud  pressure  coefficient.  It  has  been 
evaluated  by  means  of  integration  by  parts,  along  with 

Q 

the  Glauert  integrals  of  Eqs.  (2.12)  and  a  related  integral-". 
The  regular  part  may  be  expanded  cmd  written  exactly  as 
before.  Thus, 

I  A  K  d^  =  -  {0 [0„]  (a^) 

J  m  m  V  s  Km 

-X 

X 

I  dx  =  -  {0 [0„]  {b^} 

j  m  m  V  s  Km 

-X 

[Oj^]  =  [j]  -  [q]  (3.7) 


The  matrix  [j]  comes  from  the  singular  part  and  is 
given  by. 


i 


o  o 


o  o 


o  o 


c 


o 
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00 

oo 


O  O  O  /^IlTx 

O  O  o 

O  ^loo  o 

I 

xioj  o  ,<!^  o 

I 

O  r<\^>  O  O 

I 

,<1^  o  o  o 

I 

o  o  o  o 

III 

*-) 
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The  matrix  [q]  from  the  regular  part  is  analogous  to 
[p]  in  Eq.  (3 •3)  with  the  elements  ^  replaced  by 
^  in  Eq.  (2.21),  the  Fourier  coefficients  ^  by 
in  Eq.  (2.10),  and  by  in  Eq.  (2.7). 

To  complete  the  matrix  form  of  Eqs.  (1.10),  we 

set 


''rm  ^  '‘r'm  =  f^r.r') 

®r'm  “  (3-S) 

where  {a^  p/)  and  (bp/j  are  column  vectors.  The  elements 

are  derived  from  Eqs.  (2.14)  using  (A^,  +  A,,/  )  and 

I  m  I  m'  1  m 

respectively  instead  of  . 

From  these  results  and  Eqs.  (3.5)  and  (3.7).  Eqs.  (I.IO) 

give 


(bp.)  =  -  +  0  [Oj^]  (a_^)  (3.10) 

By  suitable  manipulation,  we  may  solve  and  obtain, 

^  tar.r'J  ^  ° 

■"  '  Ic^-l  +  o®  to^l 
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-[0-.]"^  {bp/}  -  a  [0  /]"^  {ap  p/} 

{b  )  =  — K - 1 - K - TjT — 

(3.11) 

if  we  assume  that  all  the  inverse  matrices  exist.  The 
matrix  [O  ]  can  be  removed  by  multiplication  with 
[O^]  ;  and  the  matrix  [O^/]”^  can  be  expanded  as  for  the 

K  lx 

zeroth  harmonic,  see  Eq.  (2.^0). 

Eqs.  (3.11),  then,  constitute  the  solution  by  this 
method.  The  homogeneous  solutions  correspond  to  Eqs.  (3.10) 
with  the  terms  {ap  p/)  and  {bp/}  suppressed,  requiring 
the  determinant  of  the  common  denominator  of  Eqs.  (3. 11) 
to  vanish. 

3.^  Comparison  of  Methods 

Which  of  these  methods  is  better  for  numerical  cal¬ 
culations  is  difficult  to  judge.  The  first  method  imposes 
an  additional  differentiation  and  iiitegration  of  the  pro¬ 
peller  contributions,  but  the  second  method  imposes  the 

determination  of  the  Fourier  coefficients  for  both  R 

m 

and  k'  .  On  the  other  hand,  the  total  number  of  operations 
m 

appear  less  for  the  first  method,  but  the  second  method  in¬ 
volves  fewer  types  of  operations. 

By  either  method,  the  complexity  of  the  calculations 
demand  every  simplification  and  utmost  care.  In  particular, 
the  limiting  case  of  infinite  advance  ratio  should  be  carefully 
checked.  For  this  case  the  equations  naturally  decouple 
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and  the  kernel  is  related  by  Eq.  (1.16)  to  Weissinger's 

function  "U  "  v;hlch  has  been  computed  for  three  harmonics, 
m 

The  corresponding  Fourier  coefficients  have  been  computed 

■3 

by  the  Bureau  Technique  Zborowski-^  over  a  range  of  values 
and  the  propeller  terms  have  been  simplified  by  explicit 
evaluation  of  the  inner  integrations  over  the  wake. 


CONCLUSIONS 


An  iteration  procedure  has  been  developed  for  the 
general  harmonic  solutions  to  the  equations  for  the  ducted 
propeller  with  finite  blade  number  in  uniform  motion  at 
zero  incidence.  Solutions  have  also  been  obtained  by 
direct  inversion. 

The  results  for  the  zeroth  harmonic  or  the  average 
shroud  load  show: 


Each  iterate  is  expressed  in  identical  form  and  is 
given  explicitly  in  terms  of  the  effective  sectional 
v-cUrtber  and  the  chord  to  diameter  ratio. 

The  iteration  procedure  reduces  to  a  matrix  form 
which  permits  rapid  numerical  calculations . 

The  nximber  of  iterations  required  for  a  prescribed 
computational  accuracy  may  be  estimated  a  priori. 

The  solution  by  iteration  converges  in  the  limit 
to  the  direct  inversion  of  the  ring  airfoil  problem 
by  J.  Weissinger. 

The  distribution  of  the  average  static  pressure 
difference  across  the  shroud  is  identical  to  the 
load  on  the  equivalent  airfoil  in  two-dimensional 
flow. 


The  results  for  the  higher  harmonics  reveal: 


The  equations  may  be  decoupled  into  the  saune  form 
as  the  equation  for  the  zeroth  harmonic  and  solved 
by  iteration. 

An  alternative  method  for  solution  by  direct,  simul¬ 
taneous  inversion  appears  to  offer  more  promise  for 
future  studies . 
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The  inherent  complexity  of  these  harmonics  requires 
careful  consideration  of  their  significance  to  design 
before  their  calculation  is  attempted. 

To  assess  the  magnitude  of  these  harmonics  and  to 
check  out  any  ccanputational  program  for  them,  the 
special  case  of  infinite  advance  ratio  is  recommended 
for  initial  evaluation. 


The  iteration  procedure  may  be  applied  to  several  other 

problems;  for  example,  the  two-dimensional  airfoil  in  a 

1 0 

bounded  stream  and  the  translational  characteristics  of 
the  GEM^^.  It  affords  not  only  important  physical  insight, 
but  also  marked  simplicity. 
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